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ABSTRACT 



Differential i^'s-band luminosity functions (LPs) are presented for a complete 
sample of 1613 nearby bright galaxies segregated by visible morphology. The 
LP for late-type spirals follows a power law that rises towards low luminosi- 
ties whereas the LPs for ellipticals, lenticulars and bulge-dominated spirals are 
peaked and decline toward both higher and lower luminosities. Each morpholog- 
ical type (E, SO, SO/a-Sab, Sb-Sbc, Sc-Scd) contributes approximately equally 
to the overall i^'s-band luminosity density of galaxies in the local universe. Type 
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averaged bulge/disk ratios are used to subtract the disk component leading to 
the prediction that the i^'g-band LF for bulges is bimodal with ellipticals domi- 
nating the high luminosity peak, comprising 60% of the bulge luminosity density 
in the local universe with the remaining 40% contributed by lenticulars and the 
bulges of spirals. Overall, bulges contribute 30% of the galaxy luminosity density 
at Kg in the local universe with spiral disks making up the remainder. If bulge 
luminosities indicate central black hole masses, then our results predict that the 
black hole mass function is also bimodal. 

Subject headings: galaxies: luminosity function, mass function — galaxies: for- 
mation — galaxies: elliptical and lenticular, cD — galaxies: spiral — galaxies: 
bulges — infrared: galaxies 



1. Introduction 



Luminosity functions (LFs) are one of the key statistical instruments widely used to 
better our understanding of galaxies, because with just a few parameters, they concisely 
describe entire populations of objects whose intrinsic properties (e.g., mass, luminosity, etc.) 
vary over orders of magnitude. In the Schechter (1976) formulation, those parameters are 
M*, the absolute magnitude corresponding to the knee of the LF, 0*, the mean galaxy space 
density, and a, the faint-end slope of the LF. The basic technique is to measure and then 
compare these parameters for different populations and thereby assess, for example, how 
galaxies evolve with redshift, or as a function of environment, or indeed any other variable 



that can be controlled for. This tool has been used, as of late, to analyse sever a. 



large 



samples of galaxies reaching out to redshifts ^ ~1 


and beyond (e.g.. 


Norbere et al. 


2OO2I; 


Blanton et al. 


2003 


Wol 


' et al. 


2003;lBell et al. 


2003; ] 


Ibert et al. 


2005 


Willmer et al. 


2006; 


Brown et al. 


2007; 


Faber et al. 


2007; 


Wake et al. 


2006; 


Reddv et al. 


2008). 



Galax y morphology is a taxonomy d evised by lHubble I (119361 ) and refined by lde Vaucouleurs 



(I1959I ) and lde Vaucouleurs et al. I (Il99ll ) that is based on the relative prominence of the stel- 
lar bulge and the degree of resolution of the spiral arms. The phy sical significance of galaxy 



morphology is that it reflects the galaxy merger history (e.g., iBalland. Silk. &: Schaeffer 
19981 ). It is conventional wisdom that the gala xy LF does depend on morphology in the 



5-band. ISandage. Binggeli fc Tammann I (119851 ) found that LFs depend on morphological 



type for galaxies in the Virgo cluster and subsequent works showed that this depe ndence 



extends to field galaxies ( 


Bineeeli, Sandaee &; Tammann 


1988 


ever. 


Efstathiou, Ellis, & Peterson ( 


1988) and 


Marzke et al. 



Lovedav et al. Ill992l) . How 



1994 119981 ) offer a contrary 
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view. More rec ently, LFs based on morphological proxiejll, such as color, also show a de- 



pendence (e g..lBlanton et al. I I2OOII: Bell et al. II2004I : iBlanton II2OO6I : IWillmer et al. II2OO6 



Faber et al. 1120071 ). at least at visible wavelengths. These studies and others indicate that 
the present-day LFs are the result of a complex evolutionary history. 

By virtue of its uniformity, reliabi lity, and full-sky coverage, the Two Micron All Sky 
Survey (2MASS, ISkrutskie et al. 1120061 ) has recently made it possible to extend LF studies 
into the infrared regime. 2MASS has been exploited in recent years to produce near-infrared 
luminosity functi ons for galaxies wi th ever greater precision, facilitated by redshifts generated 
from the SPSS (lYork et al. II2OOOI) and the 2 and 6 degree Field Galaxy Redshift Surveys 



fICole et al. I I2OOII : iKochanek et al. 1 1200 ll : iBell et al. I l2003l : Eke et al. I l2005l : IJones et al 



20061 ). The 2.16/im iT^-band (hereafter i^'-band) LF is of particular interest because of its 
relevance to understanding galaxy evolution in the context of Lambda Cold Dark Matter 
(ACDM) cosmology: at z = 0, the i^-band light traces the stel lar mass accumulated in galax- 



ies at a wavelength where i nterstellar extinc tion is minimal (jPevereux. Becklin fc Scoville 



19871 : iBell fc de Jong II2OOII : iBell et al. II2003I ). avoiding strong dependence on metallicity and 
stellar population age. For this reason, i^'-band LFs complement those obtained at visible 
wavelengths. 



Both iKochanek et al. I ( 20011) and iBell et al. I (120031 ) investigated the impact of mor- 



phology on the -ft^-b and LF. iBell et a. 



( 20031 



]) used an SDSS concentration parameter 



as a proxy for type. IKochanek et al. I (l200ll ) carried out a painstaking typing procedure in 
such a way as to permit a quantitative understanding of the uncertainties involved. Both 
studies divided the galaxy sa niples into two bro ad categories; early and late. They arrived at 
consistent total K-hand LFs (IBell et al. II2003I . Figure 9), and found that the early- and late- 
type LFs were well-described by Schechter functions that differed slightly in M* and 0*, but 
had essentially identical shapes. This is somewhat surprising given the distinctions revealed 
by LFs measured in the visible. Clearly more work remains to be done on the morphological 
type dependence of LFs in general, and for the fT-band LF in particular. Fortunately, the 
time-consuming task of assigning visible morphologies has now been completed for the vast 
majority of nearby galaxies. The principal aim of this paper is therefore to use these nearby 
galaxies to defin e new benchmark i^"-ba nd L Fs for galax i es. A s such, our study offers an 
improvement on IKochanek et al. I (1200 ll ) and iBell et al. I (120031 ) by using the most recent 
distances and finely divided morphological types. 



^Morphological proxies such as color are adopted until the difficult and time consuming task of assigning 
morphologies to the ever-growing number of cataloged galaxies can be accomplished e.g., The Galaxy Zoo 
Project, http: / / www.galaxyzoo.org/| 
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The paper is organized as follows. Section 2 describes how the nearby galaxy sample is 
selected. A non-parametric (Choloniewski) method is used to generate the LFs as described 
in section 3. The results, presented in section 4, include the morphological type dependence 
of the i^'-band LFs as well as the contribution of each to the luminosity density in the 
local universe. Recently published data have also provided the opportunity to compute the 
bulge luminosity functions and hence the contribution of bulges and disks to the i^-band 
luminosity density in the local universe. These results are also presented in section 4. A 
discussion, in section 5, explains how the results may be used to constrain models of galaxy 
evolution and includes a prediction that the bulge LP and by association, the black hole 
mass function, is bimodal and depends on Hubble type. Conclusions follow in Section 6. 



2. Sample Selection 

The goal of this study is to quantify the dependence of the K-hand LFs on the visual 
appearance of galaxies. This dictates that the sample be composed of nearby galaxies which 
have the most reliable morphological assignments. Nearby galaxies were identified using 



HYPERLEDA; a web-based interface (http://leda.univ-lyonl.fr) that provides access to the 



Principal Galaxy Catalog (hereafter PGC, Paturel et al. II2003I ). The PGC is a homogeneous 



database of galaxy parameters in the sense that an attempt has been made to place indepen- 
dent measures on a standard system. Parameters employed in this study include coordinates 
accurate to ~2", visible morphological T types and recession velocities for all known galaxies 
brighter than B = 18 mag. The PGC is a dynamic resource that is constantly being updated, 
but as of 2008 January there were 7406 nearby galaxies with Vg^^ < 3000 km s~^, apparent 
total blue magnitudes m^ < 18 mag and Galactic latitudes |6| > 10 degrees. The latitude 
constraint was imposed to avoid the inevitable incompleteness due to obscuration within our 
galaxy. However, with the exception of the Galactic plane, the PGC galaxies from which the 
sample is selected are distributed over the entire sky thereby minimizing the effect of cosmic 
variance due to large scale structures present in the sample volume. 

The 2MASS counterparts of the PGC galaxies were identified on the basis of positional 
coincidence. For an association to be made, the PGC J2000 coordinate had to fall within 10" 



of the 2MASS J2000 coordinate listed in the Extended Source Catalo£n(XSC, Jarrett et al. 



2000l ). This comparison yielded 5034 detections at 2.2 fim (i^-band) corresponding to a 



^Final checking revealed 30 galaxies with PGC and XSC coordinate differences greater than 10". The 
discrepancy is due to different centroids in visible and infrared light which is a problem particularly for 
nearby galaxies of large angular size. These galaxies are included in Table 1. 
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detection rate. Figure 1 shows the distribution of apparent i^'-band isophotal magnitudes 
for nearby galaxies (Vg^r < 3000 km s~^) in the 2MASS XSC The distribution turns over 
a.t K = 10 mag indicating that the XSC is incomplete for galaxies fainter than that. On 
this basis, a volume-limited sample was defined for further study, hereafter the KlO/3000 
sample, comprising 1613 galaxies with K < 10 mag, Vgsr < 3000kms~\ and |6| > 10 degrees 
(Table 1). The choice to use Vgs,. allows us to define a spherical volume centered on the 
Milky Way that simplifies the LF calculation. 

Every attempt has been made to identify all known galaxies with K < 10 mag, Vgsr 
< 3000 kms~^, and |6| > 10 degrees. Nevertheless, sample incompleteness can arise in two 
ways: 1) galaxies that are in the 2MASS XSC but have no radial velocities in the PGC or 
are missing from the PGC altogether, and 2) PGC galaxies missing from the XSC. A firm 
upper limit on the former can be set by noting that there are 388 objects with K<10 and 
|6| > 10 degrees in the 2MASS XSC for which no association with a PGC galaxy can be 
found. Inspection of 2MASS images reveals that most are Galactic star clusters (or indeed in 
some cases small parts of star clusters treated by the automated XSC extraction procedure 
as distinct objects), but some are galaxies. Redshifts have been determined by J. Huchra, L. 
Macri, T. Jarrett, J. Mader, A. Crook, R. Cutri, T. George, N. Martimbeau, S. Schneider, 
& M. Skrutskie (ApJ Supp., 2009, in preparation) for all 166 of the galaxies. All but three 
of which have radial velocities >3000kms~^. The three with V < 3000 kms~^ should in 
principle be included in the KlO/3000 sample but are not because they have no entries in 
the PGC; this gives an estimated completeness for Table 1 of 99.8%. With regard to the 
second source of incompleteness, there are two well known PGC galaxies of large angular 
size missing from the 2MASS XSC: the Large and Small Magell anic Clouds. These and 



some other local group dwarf galaxies are not included in the XSC (jJarrett et al. 1120031 ) and 
hence are not included in the KlO/3000 sample. However, it is unlikely that there are any 
PGC spiral and elliptical galaxies missing from the KlO/3000 sample because the = 10 
sample limit corresponds to -B = 13.5 for a normal B — K = 3.5 galaxy color, which is 4.5 
magnitudes brighter than the PGC catalog limit. Thus, only galaxies with very blue c olors 



and/or very low surface brightness would be missing from our sample (lAndreon II2002I ) and 
both selection effects will tend to further bias our sample against late-type (Sd and later) 
spirals and dwarf irregular (Im) galaxies, as explained in more detail in the next section. 



2.1. Morphological Types in the KlO/3000 Sample 

The observed distribution of K magnitudes for galaxies in the XSC with Vgsr < 3000 km s~^ 
is compared with the distribution of blue magnitudes for the entire inventory of 7406 PGC 
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galaxies contained within the same volume (Figure 1). The figure shows that 2MASS de- 
tected essentially the same number of nearby galaxies with K< 10 mag as would be found 
in a i?-band selected sample with B< 13.5 mag. 

Figure 2 illustrates how the B — K color distribution for a i^-band limited sample differs 
from a 5-band selected one. Although a 5-band selected sample with B< 13.5 mag contains 
about as many galaxies as a i^-band selected sample with < 10 mag, the latter sample 
excludes blue galaxies contained between the two magnitude limits, figuratively speaking, 
between 6 and 8 o'clock on the plot. However, the deficiency is almost exactly compensated 
for by the inclusion of additional red galaxies between 12 and 2 o'clock. Thus, a i^-band 
magnitude limited sample will contain more red and fewer blue galaxies compared to a 
5-band selected sample of similar size. 

Since a i^'-band selected sample contains fewer blue galaxies compared to a 5-band 
selected sample of comparable size, one can therefore anticipate that morphological types in 
the KlO/3000 sample will contain fewer star forming late- type spirals and dwarf irregular 
galaxies, compared to a S-band selected sample of similar size. Figure 3 supports this 
expectation. The distribution of morphological types in the KlO/3000 sample spans the 
entire range from ellipticals to dwarf irregulars. However, compared to the sample of PGC 
galaxies with B < 13.5 mag, the KlO/3000 sample appears to be missing more than 50% 
of galaxies with morphological types T > 6, corresponding to Sd and later. Consequently, 
these types of galaxies will not be considered further as there are too few to reliably define 
a LF. On the other hand, the K-hand is sensitive to the red luminous mass component in 
all galaxy types with the result that the KlO/3000 sample contains 80% of the Sc-Scd types 
and all the earlier types, plus more, than would be found in a 5-band selected sample of 
comparable size. 

A completeness test that makes no assumptions concerning the galaxy distrib ution and 



is unaffected by the presence of large scale structure was proposed by lRauzy I (1200 ll ). The test 
assumes that the LF of the population does not depend on the three-dimensional redshift- 
spatial distribution and that the apparent magnitude limit can be described by a sharp 
cutoff. In practice the test examines the distribution of the random variable 



which measures the ratio between the integrated LF up to the absolute magnitude of a 
given galaxy and the total range in magnitudes accessible at the distance of the galaxy. 
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Completeness is measured by the variable 

y ^ ESr(o-o-5) 

where the ^ represent the variance of the Q estimators. The variable has an expectation 
value of zero and unit variance. In practice the test is evaluated for sub-samples selected 
at progressively fainter apparent magnitude limits. The Tc statistic will fluctuate around 
the value of zero, but becomes systematically negative once the sample becomes incomplete. 
From the behavior of this statistic the limiting magnitude can be inferred. Figure 4 illustrates 
the result of applying this test to the full KlO/3000 sample and to several sub-samples 
segregated by morphological type. Both the number counts and the Rauzy tests suggest 
that the current sample is not affected by incompleteness for the morphological types we are 
considering, namely ellipticals, lenticulars and spirals up to and including types Scd. 



2.2. Galaxy Distances 



Because the sample galaxies are nearby, their peculiar velocities can be a significant 
fraction of their Hubble flow velocities, particularly for those with Vg^r < 1000 km s~^. Ma- 
jor perturb ers include the Virgo clus ter, which lies inside the sample volume, and the Great 



Attractor (ILynden-Bell et al. I Il988l ) which lies outside the volume but still perturbs the 



Hubble flow within. Th us, one can no t simp ly deduce a distance from the recession velocity 
and a Hubble constant. iTuUy et al. have recently quantified the peculiar velocities 

for nearby galaxies utilizi ng redshift-independen t distances based on the following m ethods; 
the Tully-Fisher relation (ITully fc Fisher 1119771 ). Ce pheids (IFreedman et al. 1120011). t he lu- 



minosity of stars at the tip of the red giant branc h (IKarachentsev et al. 



2004 



surface brightness fluctuations (ITonry et al. 11200 ll ). These methods have yielded quality dis- 



20061 ). and 



tances, with distance modulus uncertainties < 0.1 mag, for 591 nearby (Vgs^ < 3000 km s 
galaxies of which 302 are included in the KlO/3000 sample. Tully determined distances 
for other galaxies based on associations with groups which contain o ne or more n i ember s 
with quality distances together with the Numerical Action Models of IShaya et al. I (119951 ). 



Collectively, datasets provided by Tully (2007, private communication) yielded distances for 
1575 or 98% of the 1613 objects in the KlO/3000 sample and they are listed in Table l|f|. 



■^See also The Extragalactic Distance Database (EDD) at 



http : / /odd . if a. hawaii . edu/ 
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2.3. 2MASS magnitudes 



In the following analysis, isophotal magnitudes (measured within the 20 magarcsec"^ 
elliptical isophote; the parameter: k_m_k20f e in the 2MASS XSC) are adopted to charac- 
terize the dependence of f^-band LFs on galaxy morphology (section 4.2). The main reason 
for this choice is that the extrapolated total magnitudes pubhshed in the 2MASS XSC are 
unreliable for elliptical galaxies due to a restriction on th e choice of Sersic index that caused 
their total flux to be underestimated (ILauer et al.l 120071 ). by an amount we determine to be 
~0.3 maj^. Ad ditionally, extrapolated tot al magnitudes in the 2MASS XSC can lead to un- 



physical colors (jKarachentsev et al. II2002I ) suggesting that the extrapolations are unreliable 



for some spiral galaxies as well. With the exception of the ellipticals, the mean difference 
between the isophotal and total magnitudes cited in the 2MASS XSC is small: about 0.14 
mag for the bright, K < 10 mag, galaxies in our sample and is independent of galaxy type for 
morphologies spanning SO to Scd. Thus, the isophotal magnitudes are adopted as published 
with no corrections with the understanding that they slightly underestimate the total magni- 
tudes for spiral and lenticular galaxies but significantly underestimate the total magnitudes 
for ellipticals: a detail that is addressed in sections 4.3 and 4.4. 



3. Luminosity Function Determination 



The LP calculation uses the maximum likelihood method of ICholoniewski I (119861 ) which 
assumes that the spatial distribution of galaxies and their luminosities are uncorrelated. The 
merit of this common Poisson assumption is discussed further in the Appendix. By counting 
galaxies in a plane defined by distance and luminosity it is possible to obtain, simultaneously, 
the density distribution as a function of d istance and the prop e rly normalised LP, unaffected 



by de nsity variations within the sample (ICholoniewski Ill986l : iTakeuchi. Yoshikawa fc Ishii 
2000l ). The Choloniewski method is non-parametric, i.e., it makes no assumptions of a 



functional form for the LP and it yields the overall normalization (galaxy space density) 
directly. A disadvantage is that it requires binning the data. These characteri stics dis- 



and 


Efstatl 


liou. Ellis. & Peterson 


( 


1988) 


(1997h and 


Takeuchi. Yoshikawa & Ishii 



The interested reader is referred to Willmer 



''The average difference between the isophotal and total magnitudes quoted in the XSC is ^ O.f mag for 
KlO/3000 elliptical galaxies. However, based on an analysis described in section 4.4, we find that the total 
magnitude is likely to be 0.3 mag brighter, on average, than the total magnitude quoted in the XSC for 
KlO/3000 elliptical galaxies. 
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methods used to calculate galaxy LFs. Generally speaking, the results obtained using the 
Choloniewski method always agree, within the estimated stat i stical u ncertainti e s, wit h the 
other methods described by iBinggeli. Sandage fc Tammann I (119881 ): IWillmer I (119971 ) and 
Takeuchi. Yoshikawa fc Ishii I (120001 ) . The same is true for the KlO/3000 sample considered 
here. 

The Choloniewski method is implemented by plotting the two independent quantities; 
distance modulus and absolute magnitude, as illustrated in Figure 5. Galaxies are then 
binned and summed, vertically and horizontally. A small penalty is incurred as a result of 
the binning procedure because galaxies contained in partial bins, that are bisected by the 
apparent magnitude limit of the survey, must be excluded. Additionally, a few galaxies are 
excluded by the upper and lower bounds of the absolute magnitude and distance modulus 
limits. Table 2 provides a summary of the binning parameters used to generate the various 
LFs presented in this paper. The summa tions are used to itera tively solve the simultaneous 
equations 18, 19, and 20 (or 21) cited by I Choloniewski I (119861 ). These non-linear equations 
converge surprisingly quickly to yield 0(M), the differential LF, where M is absolute mag- 
nitude, p(/i), the number density of galaxies as a function of distance modulus, /i, and n, 
the average number density of galaxies in the sample. The procedure was applied to the full 
KlO/3000 sample and also to subsets of galaxies sorted by morphological type. A detailed 
explanation of how to correctly populate the covariance matrix, from which the statistical 
uncertainties are derived, is provided in the Appendix. 

In order to test the LF calculation, it was applied to samples of si ze comparable 



to th e KlO/3000 sample drawn from the Millenium Simulation database (jSpringel et al. 



Croton et al. 



2005h and in partic ular the semi-analytic galaxy catalog within it (jPe Lucia fc BlaizotI 



2007 



20061 ). For 16 independent samples drawn according to the n and rriK selection 



of the KlO/3000 sample, the derived LFs agree well with the true one derived by counting all 
simulated galaxies within the volume. The dispersion among the 16 realizations is, however, 
about 50% larger than the calculated LF uncertainties, even after normalizing the 16 sam- 
ples to eliminate cosmic variance. The excess dispersion is also about 50% larger than the 
Poisson uncertainties. Thus, to the extent that the Millenium Simulation is representative 
of the distribution of galaxies in the local universe, neither the Poisson estimate nor the 
Choloniewski estimate is a good representation of the actual uncertainties as measured by 
the standard deviation of the simulated LFs. The most likely reason for the discrepancy 
is that real galaxies are clustered. Clustering will increase the uncertainties because the 
statistically-independent unit consists of multiple galaxies (on average 2.2 of them if the 
uncertainties are increased by a factor of 1.5). The effect of clustering is discussed in more 
detail in the Appendix. For the Choloniewski method in particular, the discrepancy between 
the calculated uncertainties and their Poisson values is largest for elliptical, lenticular, and 
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early-type spirals (Table 3) which tend to be the most clustered. The Choloniewski un- 
certainty estimates thus seem particularly sensitive to clustering, but future improvements 
may be possible by implementing the method with a generalized Poisson distribution in- 
stead. Overall, the simulation results suggest that the LFs are reliable but the calculated 
uncertainties are underestimated by about a factor of 1.5 and using the Poisson method to 
estimate uncertainties would not improve the results. 



4. Results 

4.1. Density Function for the KlO/3000 Sample 

Figure 6 illustrates the solid-angle averaged radial number density of galaxies at each 
distance in the sample volume. Since the solid angle average is over very nearly the whole 
sky (S.Stt steradians) individual structures can account for only a small fraction of any 
density peak. For example, the Virgo cluster, which is the largest structure in this volume, 
constitutes only ~ 14% of the total galaxy density enhancement seen at /i + 51ogio/i ~ 
30.5 mag. That peak is likely due to the combination of the Local Supercluster (of which 
Virgo is part) and the Southern Supercluster (a Southern hemisphere structure that includes 
the groups of Dorado and the Eridanus and Fornax clusters). Consequently, one is cautioned 
against identifying density enhancements with ind ividual known galaxy clusters or groups. 



Clear evidence for a morphology-density relation (iDressler Ill980l ) is therefore absent from 
the data. Late-type spirals (Sc-Scd) have the highest average number density in the sample 
volume and ellipticals the lowest. 

The average density of galaxies appears to trend downwards at ~ 30Mpc (/i + 51ogio/i 
= 32.38 mag, h = 1), the maximum distance of the sample. The trend is not due to incom- 
pleteness for three reasons. First, the Rauzy test results described in Section 2.1. Second, 
V/Vmax = 0.7 in the outer shell (31.88 < /i + 51ogio/i < 32.3 8 mag) indicating that the sam- 



ple is complete up to the distance limit of the survey. Third, ICholoniewski I (119861 ) reported 
the same trend using an independent sample based on the CfA redshift survey. That sample 
revealed that the density increases again farther out at ~ 80Mpc {h = 1) (/i + 51ogio/i = 
34.51 mag). Thus, the downward trend at the periphery of the KlO/3000 sample is judged 
to be real and not symptomatic of a selection effect. 

The decrease in galaxy density at 30Mpc {h = 1) is large b ut within the plausible 



range of cosmic variance. An analytic estimate of cosmic variance ( iDavis fc Huchralll982l ) 
depends on the volume sampled and the galaxy two-point correlation function. The latter is 
not known directly for the i^'-selected sample, but we assume it's the same as for visually- 
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selected galaxy s amples. A correlati on function based on the fluctuation power spectrum 
from WMAP-1 (ISpereel et al.l |2003| ) was extrapolated to the present via the method of 
Seljak fc Zaldarriaga Jl996h then transformed via a spherical Bessel function and smoothed 



with a 1 Mpc radius (D. Eis enstein, private communication, 2006). The volume integral 



[e.g. 



Newman fc David 120021 . equation 1) was then evaluated via a Monte Carlo approach. 
The resulting cosmic variance uncertainty is ~12% for the full sample and ~18% for a 
sample half a magnitude less deep (Figure 6). The true correlation function is likely to be 
larger for if-selected galaxies, which are predominantly early-type, than for visually-selected 
ones, more of which are late-type (Figure 3). This suggests the preceding estimates may be 
too small. An alternate estimate for cosmic variance makes further use of the Millennium 



Simu l ation (ISpringel et al.ll2005l ) and the semi- analytic galaxy catalog (jPe Lucia &: Blaizot 



200?! : ICroton et al.lbooeh . For this calculation, simulated galaxies were simply counted in 



64 independent spheres of 30 and 23.8 Mpc {h = 1) radii, and standard deviations were 
26% and 33% respectively of the mean galaxy density. These values would make the dip in 
the final bin of Figure 6 (or the excess in pr eceding bins) only a bout 2a. The agreement of 
our overall luminosity function with that of I Jones et al. I (120061 ) (section 4.2) suggests that 
cosmic variance is not playing a large role in our results. Regardless of its magnitude, cosmic 
variance represents a single uncertainty for the entire LF, not an independent uncertainty in 
each luminosity bin, under the assumption that the LF is independent of location. 



4.2. Parametric Fits to the K-hand Luminosity Functions by Hubble Type 



Figure 7 presents K-band LFs with Schechter function fits (jSchechter 1 119761 ) for all 
galaxies and subsets segregated by morphological type. The LFs are defined in Table 3 and 
the corresponding Schechter function fit parameters are listed in Table 4. For the total K- 
band LF, value s for the parameters M^,— Slogm/t, (j)*lh? and a agree within 2a of previous 
determinations ( 



[ones et al. 2006; Eke et al 



2005: Bell et al. 2003; Kochanek et al. 2001 



Cole et al. Il200ll ). We find a ~ 40% high e r spa ce density of galaxies in the range —23 < 



Mr-— Slogmfe < —21 mag than I Jones et al. I (120061 ) which represents a difference of about 2a. 



Thus, unlike Jones et al. ( 20061 ). we do not find a residual with respect to the best fitting 
Schechter function over that magnitude interval. 

Figure 7 shows that the i^-band LF for ellipticals is represented by a Schechter function 
that declines toward both high and low luminosities, akin to the red elliptical sequence of 
(120071 ). Number counts decrease for low- luminosity elliptical galaxies until the 



Driver et al. 



low- luminosity upturn at Mr-- Slo gm/t > —21 mag which represents the onset of the dwarf 
elliptical sequence at Mb > —18 (ISandage. Binggeli &: Tammann Ill985l ). akin to the blue 
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elliptical sequence of iDriver et al. I (120071 ). 



The LFs for lenticular galaxies (SO) and bulge-dominated early-type spirals (SO/a - 
Sbc) are very similar and therefore have been combined. Like the ellipticals, the combined 
luminosity function for lenticulars and bulge-dominated spirals is represented by a Schechter 
function that declines toward both high and low luminosities. In contrast, the LF for late- 
type spirals (disk-dominated; types Sc-Scd) is completely different than found for the other 
galaxy types. Although this LF is also is well-represented by a Schechter function, the 
function is essentially a power law over the range of luminosities for which the LF can be 
defined, with a slope a =—1.4 that predicts an increasing space density of low-luminosity 
late- type spirals with no evidence of a low luminosity turnover prior to M;^— 51ogio/i < 
— 19.75 mag. An exponential turnover is expected at the high luminosity end. 

Our findings with regard to the i^'-band LF are therefore threefold. First, we find a 
total LF that is consistent with previous work. Second, each of the three galaxy classes 
(ellipticals, bulge-dominated spirals, and disk-dominated late type spirals) have a LF with a 
distinct shape, and none of them mimics the shape of the total LF. Third, ellipticals dominate 
the space density at high luminosities, a result that is only accentuated if total magnitudes 
are considered (section 4.4), whereas late- type (Sc - Scd) spirals dominate the space density 
at low luminosities. Lying between these two extremes are the lenticular galaxies and the 
bulge-dominated spirals (SO/a - Sbc). 

The only o ther study to hav e explo red the morphological type dependence of i^'-band 



LFs is that of iKochanek et al. I (120011 ). Using visual morphological classifications, they 
segregated their sample of bright, nearby galaxies {K20 < 11.25 mag, cz > 2000 km/s) into 
just two broad cate gories; early and l ate. F igure 8 illustrates the agreement between our 
results and those of iKochanek et al. I (I2OOII ) when our sample is divided the same way as 
theirs, at T = —0.5, such that the early type sub-sample includes elliptical and lenticular 
galaxies and the late- type sub-sample includes all galaxies classified SO/a and later. The 
shape of the LFs are similar though not identical when galaxies are divided into these two 
broad categories. It is not until galaxies are more finely segregated that the differences 
between the LFs for the morphological types emerge. 



4.3. fC-band Luminosity Density by Hubble type 

Calculating the luminosity density is important as it can provide a constraint on the 
mass density of stars in galaxies, given a mass to light ratio. Parameterizing the LFs allows 
the luminosity density j to be calculated by integration. 
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0(M)lO°-^(^®-^)dM 



(3) 



where i s the absolute magnitude of the Sun, corresponding to 3.32 mag at K 
( jBell et al. II2003I ). The total i^-band luminosity density was calculated by integrating equa- 
tion 2 over the interval —25 < M/^— 51ogio/i <— 19mag using the Schechter function illus- 
trated in Figure 7. This yields (5.8 ± 1.2) x 10® h Mpc A Monte Carlo method was 
used to calculate the uncertainty. O ur value for the -ftT- b and total lumino s ity density agrees 



2001 



well with previous det erminations (jJones et al. I l2006l : iBell et al. I l2003l : iKochanek et al 



Cole et al. 1120011) . all of which are about a factor of two lower than the value reported 



by iHuang et al. I (120031 ) . See iBell et al. I (120031 ) for a discussion of the various luminosity 



density estimates in the literature. 

Elliptical galaxies (Figure 7) contribute ~16 ± 3% of the total fC-band luminosity 
density of galaxies in the local universe. This value is revised upwards to ~18% if total rather 
than isophotal magnitudes are considered (section 4.4). Lenticulars and bulge-dominated 
spirals (Figure 7) contribute ~68±14% of the total, or ~22 ± 4% for each sub-group (SO, 
SO/a-Sab, Sb-Sbc). Finally, the late- type spirals (Figure 7) contribute ~16 ± 3% of the 
total. The results are summarized in Table 5. To a good approximation, each Hubble type (E, 
SO, SO/a-Sab, Sb-Sbc, Sc-Scd) contributes equally to the overall K-band luminosity density 
in the local universe. 



4.4. i^'-band Bulge Luminosity Functions by Hubble type 



In addition to the LFs for entire galaxies, the LF of bulges alone can provide information 
on how these structural components formed. For lenticular and spiral galaxies, it is possible, 
in principle, t o decompose each galaxy into its bulge and disk components by modeling imag- 



ing data (e.g..|Peng et al. II2002I : ISimard et al.ll2002l ). An alternative approach, employed by 



Graham &: Worley I (120081 ) . is to compute the bulge a r id dis k luminosities from Sersic bulge 
and exponential disk parameters. I Graham fc Worley I (120081 ) have done this for a selection of 
about 400 nearby galaxies with i^-band imaging in the published literature, many of which 
are in the KlO/3000 sample. Thus, the resulting bulge/total luminosity ratios can be used 
statistically to predict bulge LFs for different Hubble types in the KlO/3000 sample. 

The total magnitudes adopted for spiral and lenticular galaxies are those reported in the 
2MASS XSC under the parameter: k_in_ext. The total magnitudes were then corrected to 
bulge magn itudes using the bulge/tota l luminosity ratios as a function of morphological type 



reported in iGraham fc Worley I (120081 ). Galaxies exhibit a range of bulge/total luminosity 
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ratios even within a particular morphological type. Following iGraham fc Worley I (120081 ) . the 
adopted la ranges for mbuige — mtotai are 1.12 to 2.18 mag for the SO-SO/a galaxies, 0.58 to 
2.49 mag for the Sa-Sab galaxies, and 1.23 to 3.25 mag for the Sb-Sbc galaxies. A plausible 
bulge magnitude was estimated for each galaxy using a Monte Carlo method that applied 
a type specific rabuige — ^totai correction to the total magnitude of each galaxy, randomly 
selected from a uniform distributioiifl of real numbers in the ranges quoted above. A bulge 
LF was then computed using the Choloniewski method. This procedure was repeated 55 
times leading to an average LF and an associated standard deviation. The mean bulge LF 
and the la standard deviation from the mean is reported in Table 6 for each morphological 
type. 

Ellipticals do not have disks and so there is no mhuige — mtotai correction. However, the 
2MASS isophotal magnitudes for the ellipticals underestimate the total magnitudes due to 
the fact that a restriction was ir nposed on the valu e of the Sersic index when the 2MASS 
total magnitudes were computed (ILauer et al.ll2007l ). We determined the value m^so — ^totai 
= 0.44 ± 0.28 mag by comparing the 2MASS XSC i^-band is ophotal magnitu des with the 



total fT-band magnitudes measured by lMarconi fc Hunt I (l2003l ) using GALFIT (jPeng et al 



20021 ) for 14 ellipticals in the KlO/3000 sample. We found no dependence of mj^o — mtotai 
with absolute i^'-band magnitude. Thus, total magnitudes are estimated by applying a 
correction, randomly selected from a uniform distribution of numbers in the range —0.16 to 
—0.72 mag, to the isophotal magnitude of each elliptical in the KlO/3000 sample. Then a 
LF was computed using the Choloniewski method. This procedure was repeated 55 separate 
times leading to the average LF and the associated la standard deviation reported for 
ellipticals in Table 6. 

The derived bulge LFs, illustrated in Figure 9, predict that the i^-band bulge LF is 
bimodal with the ellipticals clearly offset from the bulges of lenticulars and SO/a - Sbc spiral 
galaxies. This is a consequence of the type dependent magnitude corrections employed to 
obtain the bulge LFs. The corrections essentially shift the isophotal LFs (Figure 7) towards 
higher luminosities for the ellipticals and towards lower luminosites for the lenticulars and 
the so/a - Sbc spiral galaxies, thereby accentuating the difference between the types already 
noted in section 4.2. Our combined bulge LF is similar in shape and amplitude to the B- 



band one reported by [Driver et al. I (120071 ). although a detailed comparison is complicated 



by the different wavelengths. The dwarf elliptical sequence appears to connect onto the 



^ It is not clear which distribution to adopt because indulge — '^totai has been measured for only a few 
galaxies in the K-hand. However, similar results are obtained if a normal distribution is used instead of a 
uniform one. In either case, the mean correction affects the lateral displacement of the LF along the abscissa 
and the range of the correction affects the uncertainty which is reflected in the width of the band (Figure 9). 
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LFs for spiral bulges. However, our LF for dwarf ellipticals is incomplete, these objects 
have a distinct LP that continu e s to r ise towards low luminosities as shown previously by 
Sandage. Binggeli fc Tammann I (119851 ) . 



Integrating the bulge LFs using equation 3 allows an estimate of the contribution of 
bulges to the i^-band luminosity density in the local universe. The results are presented 
in Table 7. Interestingly, bulges in each of the Hubble types, SO-SO/a, Sa-Sab, and Sb- 
Sbc contribute approximately equally to the luminosity density. Ellipticals, on the other 
hand, contribute about 50% more to luminosity density than the lenticular and spiral bulges 
combined. Collectively, all bulges contribute 30 ± 7% of the total i^'-band luminosity density 
of galaxies in the local universe, thus the remainir ig 70 ± 7% is attribut e d to d isks. A similar 
result was obtained previously in the l^-band by ISchechter fc Dresslerl (119871 ) . 



5. Discussion 

5.1. Morphological Type Dependence of the i^-band Luminosity Function 



Our principal new result is that the i^-band isophotal LFs depend significantly on galaxy 
morphology. Broadly speaking, the i^-band LFs for galaxies manifest in two varieties, both 
of which can be described mathematically by a Schechter function. The LFs for ellipticals, 
lenticulars and bulge-dominated spirals (SO/a - Sbc) are peaked with a fall-off at high and low 
luminosities. Although the functional forms for various types of bulge dominated galaxies are 
similar, the ellipticals are displaced about one magnitude brighter in luminosity compared 
to the lenticulars and bulge-dominated spirals, a difference that increases to ~1.4 mag if 
total magnitudes are considered (section 4.4). In contrast, the disk-dominated late-type 
spirals (Sc-Scd) are distributed according to a power law with lower luminosity systems 
more numerous than any other galaxy type. 



Overall, ou r resul ts agree qualitatively with lBinggeli. Sandage fc Tammann I (119881 ) and 
Loveday et al. I (119921 ). who showed that the 5-band LFs are different when galaxies are 
segregated by visible morphology. However, our results differ in detail, particularly with 
regard to the Sc galaxies. Subsequently, little evidence has been found for type dependent 
diffe rences in the B , SPSS- and 2MASS-band LFs for raorphological types earlier than Sc 



Scd (lEfstathiou. 



Nakamura et al.. 



llis. fc Peterson Ill988l : iMarzke et al. Ill994l . Il998l : iKochanek et al. 112001 



20031 ). However, those galaxy samples were divided into two or at most 



three broad morphological classes whereas our results are based on a finer morphological 
segregation. Another factor to consider is the reliability of morphological assignments for 
distant galaxies. Morphological classification becomes difficult for distant and faint galaxies. 
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even more so when the morphological assignm ent is based on photographic Schmidt plates. 



As noted previously by Ide Lapparent I (120031 ) . variations among the LFs in the literature 
depend largely on the criteria employed for typing the galaxies and the subsequent grouping 
of types. When the different morphologies are bundled together the distinctions between 
the LFs become less apparent as Figure 8 illustrates. Our results are based on well resolved 
nearby galaxies, so the classifications are more dependable, and the distinctions revealed by 
our analysis resulted from dividing the galaxy sample into several morphological bins. 

Given the difficulties with morphological classifications for distant galaxies one is com- 
pelled to critically evaluate the efficacy of surrogate but quantitativ e measures of galax y 



morphology, such as the SDSS color and light concentration indices ( Strateva et al 



both of which yi e 



Bell et al. 



200 



1 



d significant dif f erences in the resulting LFs (e.g. 



Consehce 2003; Bell et al. 2004; Willmer et al. 2006; Blanton 



Blanton et al. 



2001) 



2001 



2006 



Faber et al. II2007I ). Finding the ideal combination of parameters that can also be applied 
to galaxy simulations is one of the current observational challenges. However subjective the 
visual classifications may be, their principal merit is that they do lead to distinct functional 
forms for the A'-band LFs. Perhaps a useful compromise between the two approaches would 
be to classify galaxies using artificial neural networks (ANN), although this method too 
seems to have probler as with ANN clas s ificati ons leading to similar Schechter functions for 
all galaxy types (e.g., iBall et al. 1 12004| . 120061). The Zuric h Estim ator of Structural Types 



(ZEST) is another approach introduced by lScarlata et al. I ( 120071 ) that appears to do better 
by returning different forms for the LFs of bulge dominated and disk dominated galaxies that 
more closely resemble the distinctions revealed by our analysis. 

The differing LFs for bulge dominated and disk dominated systems suggests at least 
two quite distinct galaxy formation mechanisms are at work to produce the diversity of 
morphological types seen in the local universe. The next step is to establish what the 

s in the context of hierarchi- 
S emi-an alytic models 



formation mechanisms are, w hich will require modeling the L 



cal clustering scenarios (e.g., ICole et al. Il2000l : iBenson et al 



have revealed that a combina tion of cold gas accretion (jWeinberg et al. II2004I ) and feedback 



20031 



( lOppenheimer fc Dave II2006I ) can yield a mass function for galaxy disks that is similar to, 
but slightly flatter than, the slope of the halo mass function. Such models are also able 
to reproduce the peaked LFs observed for ellipticals and bulge-dominated spirals by incor- 



porating major mergers (IBarnes fc Hernquist 



19921 : iHopkins et al. Il2008br ). Work already 



underway shows that the morphological dichotomy revealed by the i^-band LFs can be 
understood within the context of galaxy evolution in ACDM cosmology (A. Benson & N. 
Devereux 2009, in preparation). 
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5.2. The Hubble Type Dependence of Black Hole Mass Functions 



Most astronomers now believe that supermassive black holes (BHs) reside in the bulges 
of, possibly, all galaxies. The principal evidence is a correlation between the bulge luminosity, 
meas ured in the near- infrared, and black hole mass ( Marconi fc HuntlboOS : Ferrarese fc Merritt 
2000). The importance of this correlation is that the near-infrared luminosity of stellar 



bulges, a component on which the morphological classification of galaxies is based, may 
be used as a surrogate tracer of the extragalactic black hole mass function (BHMF) which 
in tu rn constrains physica l models of black hole gr owth in the context of ACDM cosmol- 



ogv (IMarconi et al. 1 12004 : iHopkins et al. I l2008al Jbl: IShankar. Weinberg fc Miralda-Escude 



20091). 



Using the group 1 calibration of iMarconi fc Hunt I (120031 ) . the bulge LFs shown in Fig- 
ure 9 may be translated into BHMFs. When defining BHMFs this wa y, it is important 
to app reciate the following caveats. Firstly, the group 1 calibration of IMarconi fc Hunt 



(120031 ) is based primarily on elliptical galaxies, thus one has to assume that the same rela- 
tion applies to lenticulars and the bulges of Secondly, it has not yet been proven 
that BHs exist in all ellipticals, lenticulars and the bulges of all spirals, thus the derived 
space densities will be upper limits. Nevertheless, if BHs do exist in all bulges, and the 
same bulge luminosity - BH mass correlation holds for all, then our results predict that the 
BHMF is bimodal, and depends on Hubble type, with the ellipticals having a range of BH 
masses distinctly higher than the range of BH masses inhabiting lenticular and spiral galax- 
ies. Our total BHMF is similar to other total BHMFs in the published literature which have 
been derived u sing a variety of relationships between BH ra ass and bulge properties. As 
summarized by IShankar. Weinberg fc Miralda-Escude I (120091 ) . many of the existing BHMF 
determinations are based on total galaxy LFs, which include all morphological types, with 
subsequent corrections for the contribution of each morphological type to the total LF, and 
a further correction to estimate the contribution of the spheroid to the luminosity of each 
type. Ours is the first to be based on LFs explicitly segregated by morphological type. With 
the aforementioned caveats in mind, our results further predict that the r naiority of BHs 
reside in lenticulars and the bulges of spirals. In contrast, the BHMF of iGraham et al 



(120071 ). which is based on an empirical non-linear relationship between BH mass and bulge 
Sersic parameter, predicts a sharp decrease in the BHMF over that same interval (Figure 
9). Which of th ese two vastly different predictions is correct hinges on wh ether the bulge 



luminosity (e.g.. IMarconi &: Hunt II2003I ) or the bulge Sersic parameter (e.g.. iGraham et al 



^Sc and later typ es have not been includ ed in Figure 9 as central supermassive BHs have yet to be found 
in later type spirals (|Gebhardt et al. 1120011) . 
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20071 . and references therein) is the better predictor of the existence and mass of central BHs 



in lenticulars and the bulges of spiral galaxies. 

Given the aforementioned caveats, an upper limit to the local mass density of BHs, p., 
may be calculated by integrating the BHMFs using 

<j){M,)M.d{logM.) (4) 

- hower 

where 0(M,) is the differential BH mass function, mapped from the bulge LF (Figure 9), 
and M, is the BH mass. 

The result of the integral is very sensitive to the value of the upper mass limit which 
undoubtedly contributes to the factor of ~2 range in estimates for p, in the literature. 
Integrating between 6.5 < logmfM Rff) < 9-4, where the upper mass limit corresponds 



to the highest measured BH mass (IKaspi et. al 1 120071 ) yields a total BH mass density of 



(4.6 =fc 0.5) X 10^ Mf7) Mpc ^ (h = 0.7) in g ood agreement with previous determinations 



(IShankar. Weinberg fc Miralda-Escude II2009I . and references therein). However, our predic- 



tion that the BHMF is bimodal associates over half, (~ 60%), of the total black hole mass 
density with elliptical galaxies, the bulges of spirals making up the remainder, apportioned 
approximately equally between the lenticulars, the SO/a-Sab's and the Sb-Sbc's. 

At present, stellar bulges provide an efficient but indirect means of constraining super- 
massive black hole demographics without recourse to the difficult and time-consuming mea- 
surement of black hole masses directly using gas and star kinematics. In particular, because 
stellar bulge luminosities can be measured consistently for very large numbers of galaxies, 
the method has the potential to reduce biases that may aff ect other indirect methods, such 



as the use of broad emission lines (e.g. iGreene &: Ho 1120071 ). We are currently undertaking 



the next logical step to refine this analysis which is to isolate the bulges explicitly using the 
2MASS images. 



6. Conclusions 

i^-band LFs have been presented for a complete volume limited sample of 1613 galaxies 
with Vgsr < 3000 km s~^, K < 10 mag and \b\ > 10 degrees. Our principal conclusions are 



1. The i^-band LFs depend on morph ological type when gal a xies a re segregated accord- 
ing to the visual classification scheme of Ide Vaucouleurs et al. I (Il99ll ). Ellipticals dominate 
the space density at high luminosities whereas late-type spirals dominate the space density at 
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low luminosities. Lying in between the two extremes are the lenticulars and bulge-dominated 
spirals. 

2. The K-hand LF for late- type spirals follows a power law that rises towards low lumi- 
nosities whereas the K-hand LFs for ellipticals, lenticulars and bulge-dominated spirals are 
peaked with a falloff at both high and low luminosities. 

3. Each morphological type (E, SO, SO/a-Sab, Sb-Sbc, Sc-Scd) contributes approximately 
equally to the overall i^'-band luminosity density in the local universe. 

4. The ii'-band LF of bulges is bimodal and depends on morphological type. Ellipticals 
comprise 60% of the bulge luminosity density in the local universe, the remaining 40% is 
associated with lenticulars and the bulges of spirals. Overall, bulges contribute ~ 30% of 
the total galaxy luminosity density at K, with disks making up the difference. 

5. If the bulge luminosity traces the black hole mass, our results predict that the black 
hole mass function is bimodal and depends on morphological type with the most massive 
black holes occurring in ellipticals and the space density of black holes reaching a maximum 
in the bulges of spirals. 
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A. Appendix 



This Appendix describes in detail how to calculate the size of the error bars on the 
differential LF 0(M), the number density of galaxies p{fi), and n, the average number density, 
where M is absolute rnagnitu de and fi is the distance modulus. The reader is encouraged to 
consult ICholoniewski I (119861 ) prior to implementing the following strategy which requires a 
familiarity with the maximum-likelihood method of statistical analysis. 

In order to characterize the statistical uncertainties on 0(M), p(/i), and n, the observ- 
ables, M and /i, must be represented by a particular statistical model. Generally, this model 
is assumed to be a Poisson distribution, which means that the galaxies are distributed inde- 
pendently of position, that the average number of galaxies in a region of space is proportional 
only to the volume of that region, and that their luminosities are uncorrelated with their 
locations. This model yields the familiar and widely adopted result 1 / (2.3^/N) for the uncer- 
tainty in logio0(M), where N is the number of galaxies in the absolute magnitude interval, 
AM. However, it is rarely, if ever, proven that the observables, M and yU, are distributed 
in a Poisson fashion even though there are reasons to believe that this may not be a good 
assumption given the clustering tendencies for elliptical galaxies in particular. We have dis- 
covered, in the course of writing this paper, that some uncertainties calculated using the 
procedure described below are less than 1/2.3VN. These cases are indicated by asterisks in 
Table 3. This outcome suggests that the galaxy distribution is not Poisson, but is instead 
slightly clum ped in space and luminosity and better represented by a generalized Poisson 
distribution (IConsul Ill989l ). The generalized Poisson distribution is described by 

N{l-b). 



p{N\V,b) 



■[N{l-b)+Nb] 



N-l -N(l-b)-Nb 



(Al) 



Here p{N\V, b) is the probability that a cell of volume V placed randomly in space contains 
exactly N galaxies. If n is the average density of gala xies then N = nV. An application of this 



distrib ution to hierarchical clustering is described by lSheth I (119981 ) and lSaslaw &: Hamilton 



( ll984j ). A formal analysis, similar to that outlined below, but beyond the scope of the present 
paper, would constrain the parameter b to the range < b < 1 with the standard Poisson 
distribution corresponding to the case 6 = 0. Such a distribution would yield uncertainties 
that are larger than 1/2. by a factor 1/(1 — b). 



Following ICholoniewski I (119861) . we pr oceed with the Poisson distribution function, and 
adopt the notation of ICholoniewski I (119861 ). which is most easily understood with reference 
to Figure 2 of that paper. A is the total number of absolute magnitude columns, indexed by 
the counter i, and B is the total number of distance modulus rows, indexed by the counter 
j. A magnitude limited sample fills the M, fi plane only partially, in the region below the 
selection line M + fi = m^j^, where m^j^ is the magnitude limit of the survey. The parameter 
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S, is the magnitude limit of the survey, mum, in units of i and j, such that z + j < S, an 
inequahty because the galaxies are also contained within a finite volume and within a finite 
range of absolute magnitudes. 

The covariance, a matrix whose leading diagonal contains the uncertainties for the 
values, (f)i,pj,n, can be approximated by use of the information matrix, specifically, 



Cov{E) 



dnn{l) 



E=E 



where / is the hkelihood func t ion, w hich is chosen to be Poisson 

Jl986h 



Choloniewski 



see 



EadieJ 119821). Following ICholoniewski I (119861 ) 



(A2) 

The minus sign cor - 



rects a typo in 

E = {(f)i,pj,n), are called the estimators and are the computed values defining 0(M), the 
differential LF, where M is absolute magnitude, p{n), the number density of galaxies as a 
function of distance modulus, p, and n, the average number density of galaxies in the sample. 



The Hessian 
itself a matrix: 



of / can be written as a 3 by 3 symmetric matrix with each entry 





r d'^lnH) 






dHnil) 


a'^ln{l) 


d4>dp 
aHn(l) 


d(f>dn 
dHn(l) 


dE^ 


dpd(p 


dp^ 


dpdn 


d'^ln(l) 


d'^ln{l) 


d'^ln(l) 




dnd(f> 


dndp 





(A3) 



So, the next step in the procedure is clear, but complicated, because it entails taking the 
second derivative of the likelihood function with respect to each of the variables, (pi-, pj and 
n. 



The Poisson likelihood function is given by 



I 



(A4) 



i+j<S 

where Xij = j^(f)iAMpj{Vj — V^-i), and V is the sample volume and AM is the magnitude 
interval. 



Taking the natural log we obtain 



ln{l) 



i+j<S 



(A5) 



Unless otherwise stated double products such as ^^^-^^^j-^ will be denoted by Hi+ixs 

i+j<S ' - 

and double sums such as will be denoted by X]i+7<s- 

i+j<S 
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Here Ni^j — the number of galaxies in bin i,j. Consequently 



dcpk dpk dn 



(A6) 



For the first sub-matrix, 



dHnjl) 



d'^ln(l) 



d(l)id(j)j 



(A7) 



We have 



dnn{i) _ d 

d(f)id(f)j d(j)id(l) 



(A8) 



^ \a+3<S 



By definition Aq,/3 = ^(pa^Mpf^iVp — V^/3-1), where V is the sample volume and AM is the 
magnitude interval. We introduce new dummy variables a,j3 which are independent of the 
variables This is necessitated by the need to distinguish between the variables which 
the hkelihood is a function of and the factors which the likelihood is composed of. 

Letting = AM{Vp — Vg-i), this becomes X^p — Cp^cpaPf}, and so 



So 



d(t>id(t>j 



dXais _ if a 7^ i, 

dcpi Cf^^pfs if a = i. 

for all a, (3, Also = for all a, (3, Consequently 



Now 



d(t>id<t>j 

dHn{l) 



^ Nc,pln{Xap) . 



\a+f3<S 



d 



d(f) 



af)- 



^ \a+P<S 



a+l3<S 



a+l3<S 
dln{Cp^(paP(3) 



d(t)j 



dln{Cp'-ct>,pp) ^ 1 



13=1 /3=1 



(A9) 



(AlO) 

(All) 
(A12) 

(A13) 
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So 



Replacing the dummy variable (3 by the variable j we obtain 

d(f)id(f)j if « 7^ J 

Therefore the first block of the Hessian matrix is an A x A matrix 



(A14) 
(A15) 

(A16) 
(A17) 
















(A18) 



in which all terms are zero, except the leading diagonal, which is just the sum of the galaxies 
in each column, i, divided by -cff. 



We next calculate -q^§^- Since = Cp^(f)aPfj, 

dXaf3 ^ if (3y^J, 

dpj Cj^(f)a if P^j- 



(A19) 



52 A, 



Therefore g^.g^. — for all i,j,a, (5 and 

dHn(l) V- N 

dpidpj dpidpj ^ ^ 



(A20) 



a+f3<S 



One has 



d d ^ ^ d 1 

XI ^«/3^n(Aa^)= X N^p—ln{\a(i) = ^N^j—ln{C0-^aPi) (A21) 



a+/3<5 



a=l 



5-i 



a=l 



S-3 



(A22) 
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Therefore 



dpidpj dpi ^ pj 



aj 



ifi^j, 



if i ^ J, 
if i = j. 



(A23) 
(A24) 
(A25) 



To summarize, the next block of the Hessian matrix is a S x S matrix 



dpdp 




















(A26) 



in which aU terms are zero, except the leading diagonal, which is just the sum of the galaxies 
in each row, j, divided by -p^. 



To calculate g^i ' we proceed as follows: Since A^^ = C/s^(j)aPp, 



hence 



and 



Also 



dn 



-Cpn (paPp, 



^ = 2C,n-^<p^p,, 



- a+l3<S 



n 



a+P<S 



{ln{Cfi<\)a.p(i) - ln{n)) 



{ln{n)) = n ^. 



Consequently 



dHnjl) 



-2 ( Cp(t>^pp + I Y ^«/? " ^ 

\a-\-l3<S ) \a+^<5 



(A27) 
(A28) 
(A29) 

(A30) 

(A31) 
(A32) 



a single scalar term. 



We next calculate the first of the cross terms, ^q^q^ ■ 

Since Aa/3 = Cp^(j)aP)3, 

dXap ^ if/?7^J, 



then 



a0_ ^0 if /3 7^ i, 

ii (3 ^ j or a ^ i, 
Cj^ li (5 — j and a — i. 



Since 

dHn{Kp) 



dcpidpj 
Consequently |^ 



^ {InC, - ln{n) + ln<t>^ + Inp,) = A 



0-0 + 



dpj 



-E 



a+l3<S d<pidpj 



C. 



^ ifi + j<S, 
otherwise. 



Therefore: 



dHn{l) 
d(pdp 



-Ci/n 


-C-i/n 


-Ci/n 








-C2/n 


-C2/n 


-C2/n 




C2/n 


din 




-C^/n 






din 


din 


-d/n 


-CJn 


-d/n 




CJn 





-C^/n 













































Cb/u 

















a B X A array, and its reflection, ^q^q^ , an A x B array. 
To calculate the next cross term • 



so 



dn 



d'^Xa/3 

d(f)idn 



-Cpn (f)aPl3, 



if q; 7^ i. 
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Therefore 



d(j)idn a+i3<s ^ a^i^n 

(where x f\y denotes min{x, y)). 



= - + + Inpp) = 0- £-- + + = 0. (A42) 

So 



d(f)idn 

and 

a column vector, and its reflection, a row vector, each containing A terms. 

Now, the final cross term, 

d^ln{l) 



dpjdn ^ if/? = j. 



So 



(A43) 



(A44) 



(A45) 



dpjdn 
Since 

^ = -C^n-2(/.„p^, (A46) 
d^X^fs if/9^j, 



(A47) 



— — -Uj — -Uj 



9pj9n a+/3<5 "-^ "-^ 

(A48) 



To calculate 



dpjdn 

dlnXap d ^ , , X , , , X 1 



{InCp - ln{n) + /n^^ + Inpp) = --. (A49) 
d'^ln{Xai3) dl/n 



dn dn n 

So 

0. (A50) 



dpjdn dpj 
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Thus, 



d^ln{l) 
dpjdn 

So, the last block of the Hessian matrix is, 



(A51) 




i=l 



Aa{S-B) 




(A52) 



another column vector, and its reflection 



d'^lnjl) 
' dndp ' 



a row vector, each containing B terms. 



Finally, take the inverse of the negative of the {A + B + 1) x {A + B + 1) Hessian 
matrix, the leading diagonal of which contains the uncertainties on 0(M), and n. 
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Fig. 1. — Histograms showing the distribution of B —3.5 magnitudes for a sample of 7406 
nearby PGC galaxies (Vg^^ < 3000 km s~^, |6| > 10°, m^ < 18 mag) and the distribution of 
ii"-band magnitudes for the subset of 5034 in the same volume that appear in the 2MASS 
XSC, and from which the KlO/3000 galaxy sample is selected. 
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Fig. 2. — B— K color versus apparent B magnitude for 5034 nearby PGC galaxies in the 
2MASS XSC {Y^,r < 3000kms-^ |6| > 10", < 18 mag). Black dots highlight 1596 
galaxies in the 2MASS XSC with K < 10 mag that are contained within the same volume. 
A dashed line identifies the nominal 2MASS XSC limit at K = 14.5 mag. 
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Fig. 3. — Histograms showing the distribution of morphological types for a 5-band selected 
sample of 1711 nearby PGC galaxies (Vg^^ < 3000 km s~^, \b\ > 10°, m^ < 13.5 mag) versus 
the distribution of morphological types for 1610 galaxies in the KlO/3000 sample, contained 
within the same volume. 



- 35 - 



1 1 

_ - 1 - - - - 




' -^""""'^ 




p ^ 

\\^ 




\\ ^ 


Total 

E 

SO 

SO/a-Sab 

- - - Sb-Sbc 

Sc-Scd 
1 1 


— V\v- 
\ \. ^ 

\ N. ^ 



8 9 10 

K (mag) 



Fig. 4. — Rauzy test for 1601 galaxies in the KlO/3000 sample. The test is applied to 
subsamples of the 2MASS XSC defined at progressively fainter apparent magnitudes. In 
the present figure these levels were defined at intervals of 0.1 magnitudes. The horizontal 
dashed lines show the levels from where the null hypothesis that the sample is complete 
can be rejected at 99.9% level, under the assumption of Gaussian statistics. The solid line 
indicates the expectation value of T^. The vertical dashed line indicates the adopted limit of 
^^=10. The colored lines indicate the results for the Rauzy test for different morphological 
subsamples indicated in the figure legend. The null hypothesis can be rejected with high 
confidence beyond K ~ 10.1 where the curve for becomes systematically negative. Both 
the number counts and the Rauzy tests suggest that the current sample is not affected by 
incompleteness for the morphological types we are considering. 
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Fig. 5. — Absolute K magnitude, Mi^ — 51ogio/i, versus distance modulus, hk + 51ogio/i, 
for 1349 galaxies (black dots). Gray dots identify 226 galaxies that are excluded by the 
binning process (see text for details). The vertical and horizontal lines identify the absolute 
magnitude and distance modulus limits imposed on the sample prior to computing the total 
LF (Table 2). 
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Fig. 6. — Number density as a function of distance modulus for 1349 galaxies in the 
KlO/3000 sample segregated by morphological type. For reasons discussed in section 3 
all plotted error bars are 50% larger than calculated using the maximum likelihood method 
outlined in the Appendix. The horizontal bar under the lowest curve identifies the extent of 
the Virgo cluster which constitutes only ~14% of the total density in that distance range. 
Thus, one is cautioned against identifying density enhancements with individual structures 
because the density is angle averaged over essentially the whole sky. 
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Fig. 7. — i^-band isophotal luminosity functions for the KlO/3000 sample. Filled squares 
represent the results of this paper and solid lines represent Schechter fits to the data (squares) 
Top left panel: Entire sample of 1349 galaxies. Increasing incompleteness is expected for 
Mk — 51ogio/i > —19.75 mag due to the omission of late-type (Sd a nd later) spiral and 



dwarf irregular (Im) galaxies. Stars identify the A'-band SWML LF of I Jones et al. I (120061 ) 



for a different galaxy sample based on 2MASS isophotal magnitudes and a surface brightness 
dependent correction factor added to emulate total magnitudes. The Schecter fit excludes 
the first data point at Mk — 51ogio/i = —18.75. Top right panel: LF for 142 elliptical 
galaxies. The fit excludes the first data point at — 51ogio/i = —21.25 which is attributed 
to the dwarf elliptical sequence. Lower left panel: Combined i^'-band isophotal LF for 904 
lenticular (SO) and bulge-dominated spiral (SO/a - Sbc) galaxies. Lower right panel: LF 
for 247 late-type (Sc-Scd) spiral galaxies. The fit excludes the first data point at Mk — 
51ogio/i = —19.75. For reasons discussed in section 3 all plotted error bars are 50% larger 
than calculated using the maximum likelihood method outlined in the Appendix. 
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Fig. 8. — Comparison illus trating; the agreement be tween our results and the cz > 2000 km/ s, 
K < 11.25 mag, sample of iKochanek et al. I (120011 ) when both are divided at T = —0.5 such 
that the early-type subsample includes elliptical and lenticul ar galaxies and the late -type 
subsample includes all galaxies classified SO/a and later. The iKochanek et al. I (120011 ) LFs 
are based on 2MASS isophotal magnitudes. Plotted error bars on the squares are 50% larger 
than calculated using the maximum likelihood method outlined in the Appendix. 
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Fig. 9. — Bulge luminosity functions for galaxies in the KlO/3000 sample segregated by 
morphology. The width of each band reflects the uncertainty associated with correcting 
total magnitudes to bulge magnitudes for the SO - Sbc galaxies and correcting isophotal 
magnitudes to total magnitudes for the ellipticals. The total bulge LF merges with the 
elliptical LF near Mk — 51ogio/i = —23.3 mag. The upper abscissa and the right hand 
or dinate show the bu l ge LF s translated to black hole mass functions using the linear relation 
of iMarconi fc Hunt I J2003h : logw{MBH) = -3.757+1.13[(M0-(M;^ - 51ogio/i))/2.5] and 
logio0(MBH) = (2.5/i^/1.13)logio0(M). Black hole masses and mass functions are reported 
in solar units adopting h = 0.7 and = 3.32 mag. Dashed lines show a rang e of BHMFs in 



the literature as summarized by lShankar. Weinberg fc Miralda-Escude I (120091) . Solid squares 



show the All galaxies black hole mass function from iGraham et al. 



fl2007f ). 



Table 1. KlO/3000 Galaxy Sample 



PGC 


NAME 


2MASS ID 


RA^' 


DEC 


Kb 




Rd 


LGA° 


Morphological' 


Distance^ 


Quality'' 








Dec. Deg 


Dec. Deg. 


mag 


mag 


mag 




Ttype 


Mpc 


Distance 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


12651 


NGC1316 


03224178-3712295 


50.67 


-37.21 


5.69 


-25.97 


9.44 


Z 


-1.7 


21.48 


q 


41220 


NGC4472 


12294679+0800014 


187.44 


8.00 


5.51 


-25.66 


9.27 


Z 


-4.8 


17.14 


q 


30019 


NGC3147 


10165363+7324023 


154.22 


73.40 


7.50 


-25.47 


11.26 




3.9 


39.26 




43296 


NGC4696 


12484927-4118399 


192.21 


-41.31 


7.30 


-25.45 


11.65 


z 


-3.7 


35.48 


q 


13505 


NGC1407 


03401190-1834493 


55.05 


-18.58 


6.86 


-25.44 


10.70 


z 


-4.5 


28.84 


q 


70090 


IC1459 


22571068-3627449 


344.29 


-36.46 


6.93 


-25.40 


10.95 


z 


-4.7 


29.24 


q 



Note. — Table [T] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown hero for guidance regarding 
its form and content. The galaxies are ordered by absolute K magnitude (column 7). 

''J2000 coordinates taken from the 2MASS XSC based on the peak pixel l|jarrett et al. IboOOh . 

''K magnitude measured within the 20 mag/arc sec^ elliptical isophote, parameter k-m-k20f e in the 2MASS XSC. The photometric uncertainty in 
the magnitude is typically 3%. 

■^Absolute K magnitude computed using the distance in column 11 [h = 0.75). 

''Total R magnitude taken from the Principal Galaxy Catalog dPaturel et al. iboOSh . 

Galaxy also listed in the Large Galaxy Atlas lljarrett et al. Il2003h based on the parameter cc_flg = Z in the 2MASS XSC. 

Morphological T type taken from the Principal Galaxy Catalog jPaturel et al. II2003I) . 
^Distances provided by TuUy (2007, private communication), 
''"q" indicates quality distance (see text for details). 
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Table 2. i^-band Luminosity Function Binning Parameters 



Sample 


Mo- Slogio/J 


Ma - 51ogioA 


/io + 51ogio/i 


A 


B 


S 




mag 


mag 


mag 


mag 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Total 


-25 


-18.5 


26.5 


13 


12 


17 


Elliptical 


-25 


-21 


29.5 


8 


6 


11 


SO 


-24.5 


-20.5 


29 


8 


7 


11 


so/a - Sab 


-24.5 


-21 


29 


7 


7 


11 


Sb-Sbc 


-24.5 


-20.5 


29 


8 


7 


11 


Sc-Scd 


-24 


-19.5 


28.5 


9 


8 


11 



Note. — Following the notation of lCholoniewskil lll986l) , the apparent magnitude limit, 
^lim = 10 mag, the bin size A = 0.5 mag, the solid angle Q = S.Stt steradians and the 
upper bound on the distance modulus, fiB + 51ogio/i = 32.5, for all samples 
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Table 3. X-band Luminosity Functions 



Sample 


Mk - 51ogio/» 


\ogio<p/h^ 


N 


Sample 


Mk - 51ogio/i 


logio</>//i^ 


N 




mag 


GalMpc~3 mag~i 


Gal 




mag 


GalMpc"^ mag~^ 


Gal 


(1) 


(2) 


(3) 


(4) 


(1) 


(2) 


(3) 


(4) 


Total 


-24.75 


-3.58 ± 0.12 


12 


Elliptical 


-24.75 


-3.76 ± 0.13 


8 




-24.25 


-2.94 ± 0.06 


53 


T < -3.6 


-24.25 


-3.36 ± 0.06* 


20 




-23.75 


-2.59 ± 0.04 


120 




-23.75 


-3.20 ± 0.02* 


29 




-23.25 


-2.28 ± 0.02 


242 




-23.25 


-3.17 ± 0.03* 


31 




-22.75 


-2.18 ± 0.02 


303 




-22.75 


-3.28 ± 0.04* 


24 




-22.25 


-2.08 ± 0.02 


273 




-22.25 


-3.40 ± 0.09 


14 




-21.75 


-2.05 ± 0.03 


174 




-21.75 


-3.41 ± 0.14 


8 




-21.25 


-2.13 ± 0.04 


97 




-21.25 


-3.11 ± 0.16 


8 




-20.75 


-2.14 ± 0.08 


33 












-20.25 


-2.00 ± 0.11 


18 












-19.75 


-1.90 ± 0.14 


12 












-19.25 


-2.17 ± 0.26 


3 












-18.75 


-1.70 ± 0.22 


5 
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Table 3 — Continued 



Sample 


Mk - 51ogio/i 


GalMpc^'' mag^^ 


N 


Sample 


Mk - 51ogiofe 


logio</)//i^ 


N 




mag 


Gal 




mag 


GalMpc^'' mag^^ 


Gal 


(1) 


(2) 


(3) 


(4) 


(1) 


(2) 


(3) 


(4) 


Lenticular 


-24.25 


-3.55 ± 0.11 


13 


Sa-Sab 


-24.25 


-4.06 ± 0.21 


4 


-3.5 < T < -0.6 


-23.75 


-3.41 ± 0.09 


18 


-0.5 < T < 2.4 


-23.75 


-3.26 ± 0.07 


25 




-23.25 


-2.84 ± 0.03* 


66 




-23.25 


-2.97 ± 0.04* 


49 




-22.75 


-2.75 ± 0.02* 


81 




-22.75 


-2.80 ± 0.02* 


73 




-22.25 


-2.73 ± 0.03* 


64 




-22.25 


-2.76 ± 0.03* 


57 




-21.75 


-2.60 ± 0.05 


48 




-21.75 


-2.69 ± 0.05 


38 




-21.25 


-2.76 ± 0.09 


22 




-21.25 


-3.00 ± 0.11 


13 




-20.75 


-3.10 ± 0.22 


4 




-20.75 


-3.08 ± 0.26 


3 


Sb-Sbc 


-24.25 


-3.43 ± 0.10 


17 


Sc-Scd 


-23.75 


-3.52 ± 0.11 


14 


2.5 < T < 4.4 


-23.75 


-3.14 ± 0.06 


33 


4.5 < T < 6.5 


-23.25 


-3.16 ± 0.07 


32 




-23.25 


-2.90 ± 0.04 


58 




-22.75 


-3.01 ± 0.06 


45 




-22.75 


-2.79 ± 0.03* 


74 




-22.25 


-2.74 ± 0.04* 


55 




-22.25 


-2.60 ± 0.02* 


77 




-21.75 


-2.64 ± 0.05 


42 




-21.75 


-2.75 ± 0.06 


34 




-21.25 


-2.50 ± 0.07 


35 




-21.25 


-2.90 ± 0.10 


17 




-20.75 


-2.43 ± 0.14 


12 




-20.75 


-2.91 ± 0.20 


5 




-20.25 
-19.75 


-2.30 ± 0.23 
-2.01 ± 0.23 


5 
6 



Note. — The sum of the number of galaxies within each Hubble type does not add up to the total number of galaxies because of the 
binning procedure inherent to the Choloniewski method and the fact that every plotted point is defined by a bin containing at least 3 
objects. The uncertainties on the LF values are calculated using the maximum likelihood method outlined in the Appendix. Asterisks 
identify uncertainties that are less than l/(2.3\/lV) suggesting that the galaxy distribution is not Poisson but rather clumped in space 
and luminosity (section 3; Appendix). Galaxies ar e segre gated according to morphological T types; a numerical system developed by 
Ide Vaucouleurs ] lll959l ') and Ide Vaucouleurs et al. I lll99lf) . The low luminosity end of the total LF, — 51ogio/j > —19.75 mag, is 
defined by very late-types (T > 6.5) of which there are too few to define a LF separately. 
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Table 4. i^T-band Luminosity Function Fit Parameters 



Sample 


4>»lh^ 








M»— Slogioft 


a 




galaxies Mpc~^ 


mag" 


-1 


mag 




(1) 


(2) 








(3) 


(4) 


Total 


( 11.5 ± 3.4 ) 


X 


IQ- 


-3 


-23.41 ± 0.46 


-0.94 ± 0.10 


Elliptical 


( 17.6 ± 0.9 ) 


X 


10" 


-4 


-23.42 ± 0.17 


-0.03 ± 0.16 


SO - Sbc 


( 15.7 ± 1.4 ) 


X 


10" 


-3 


-22.49 ± 0.20 


-0.18 ± 0.16 


Sc-Scd 


( 15.9 ± 4.8 ) 


X 


10- 


-4 


-23.33 ± 0.33 


-1.41 ± 0.06 



Note. — Columns (2) - (4) refer to Schechter function parameters (Fig- 
ure 7). The uncertainties reflect unweighted fits to the data points. 
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Table 5. X-band Luminosity Density by Hubble Type 



Sample 


IO'^/iLq Mpc-3 


(1) 


(2) 


Total 


58 ± 12 


Elliptical 


8.1 ± 1.7 


SO - Sbc 


34 ± 7 


Sc-Scd 


8.0 ± 1.6 



Note. — The luminosity den- 
sties result from integration of 
the Schechter functions, defined 
in Table 4, between -25 < Mk 
< -19. 



-47- 



Table 6. Bulge Luminosity Functions 



Sample Mx — 51ogio/i 
mag 

(1) (2) 



logio(j>/h^ Sample 
GalMpc"^ mag~i 

(3) (1) 



Mk - 51ogio/i 
mag 
(2) 



logi()<^//r^ 
(3) 



Total 



-25.25 


-3.76 ± 0.09 


Elliptical 


-25.25 


-3.76 ± 0.09 


-24.75 


-3.38 ± 0.05 


T<-3.6 


-24.75 


-3.38 ± 0.06 


-24.25 


-3.26 ± 0.05 




-24.25 


-3.25 ± 0.05 


-23.75 


-3.17 ± 0.05 




-23.75 


-3.17 ± 0.05 


-23.25 


-3.31 ± 0.06 




-23.25 


-3.31 ± 0.06 


-22.75 


-3.03 ± 0.09 




-22.75 


-3.42 ± 0.06 


-22.25 


-2.73 ± 0.07 




-22.25 


-3.50 ± 0.08 


-21.75 


-2.47 ± 0.06 




-21.75 


-3.37 ± 0.08 


-21.25 


-2.39 ± 0.04 








-20.75 


-2.33 ± 0.05 








-20.25 


-2.29 ± 0.06 








-19.75 


-2.37 ± 0.08 








-19.25 


-2.69 ± 0.11 








-18.75 


-3.09 ± 0.20 
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Table 6 — Continued 



Sample 


Mk — 51ogio/i 


logio</>/h^ 


Sample M/f — 51ogioft 


logio^/ft* 




mag 


GalMpc"* mag~^ 


mag 


GalMpc"* mag~^ 


(1) 


(2) 


(3) 


(1) (2) 


(3) 


Lenticular 


-22.75 


-3.64 ± 0.09 


Sa-Sab -22.75 


-3.50 ± 0.10 


-3.5 < T < -0.6 


-22.25 


-3.27 ± 0.05 


-0.5 < T < 2.4 -22.25 


-3.20 ± 0.06 




-21.75 


-2.95 ± 0.04 


-21.75 


-2.96 ± 0.05 




-21.25 


-2.76 ± 0.04 


-21.25 


-2.90 ± 0.05 




-20.75 


-2.72 ± 0.04 


-20.75 


-2.86 ± 0.06 




-20.25 


-2.64 ± 0.05 


-20.25 


-2.88 ± 0.09 




-19.75 


-2.74 ± 0.07 


-19.75 


-3.00 ± 0.10 




-19.25 


-3.11 ± 0.14 






Sb - Sbc 


-22.25 


-3.44 ± 0.10 






2.5 < T < 4.4 


-21.75 


-3.15 ± 0.07 








-21.25 


-2.97 ± 0.05 








-20.75 


-2.84 ± 0.05 








-20.25 


-2.82 ± 0.05 








-19.75 


-2.84 ± 0.08 








-19.25 


-2.91 ± 0.10 








-18.75 


-3.09 ± 0.22 







Note. — The uncertainties reflect the ± la range in the bulge/total luminosity ratio propagated through the LFs 55 
times using a Monte-Carlo method (section 4.4). 
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Table 7. iC-band Bulge Luminosity Density by Hubble Type 



Sample 


lO^'ftL© Mpc-3 


(1) 


(2) 


Elliptical 


10.7 ± 1.3 


SO - so/a 


2.6 ± 0.4 


Sa - Sab 


2.7 ± 0.6 


Sb - Sbc 


1.8 ± 0.3 


Total 


17.8 ± 1.5 



Note. — The luminosity den- 
sties result from integration of 
the bulge LFs, shown in Fig- 
ure 9, between —25 < < 
-19. 



